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Abstract 

We discuss a formulation of the fusion procedure for integrable models 
which is suitable for application to non-standard R-matrices. It allows 
for construction of bound state R-matrices for AdS/CFT worldsheet 
scattering or equivalently for the one-dimensional Hubbard model. We 
also discuss some peculiar cases that arise in these models. 






1 Introduction 


Integrable systems constitute a special class of physical models that are exactly solvable 
Q. The key ingredient that allows for the explicit construction of exact solutions is the 
so-called R-matrix. For most known integrable models, the corresponding R-matrices 
are determined by the underlying symmetry algebra of the system. This is usually an 
inhnite-dimensional Hopf algebra of Yangian or quantum affine type. 

Computing the R-matrix in various representations of such an algebra then describes 
different types of particles. For example, the Heisenberg XXX model has the Yangian of 
su(2) as its symmetry algebra. The R-matrix in the fundamental representation simply 
describes a chain of spin-^ particles. Similarly, particles of higher spin can be considered 
by taking higher-dimensional representations of 5u(2). 

On the other hand, from basic representation theory it is well known that for example 
the tensor product of two spin-^ particles splits into a spin-1 and a spin-0 representation. 
In other words, one should be able to relate R-matrices in higher dimensional repre¬ 
sentations to the fundamental R-matrix. For example, the R-matrix of spin-1 particles 
should allow for some decomposition into the R-matrix of spin-| particles. Similarly, from 
the fundamental representation it should be possible to construct new R-matrices corre¬ 
sponding to other representations. This construction goes under the name of fusion |^. 
For most integrable systems this procedure is well-understood, however, the established 
formulas do not directly apply for some of the more exotic integrable models. 

Recently, there was renewed interest in the held of integrable systems due to their 
appearance in string and gauge theory via the AdS/CFT correspondence ( see |[^ and 
references therein). In particular, this integrable structure gave rise to an unusual R- 
matrix that displays Yangian symmetry corresponding to the centrally extended 5u(2|2) 
Lie superalgebra (^5|. Remarkably, this R-matrix turned out to be directly related to 
Shastry’s R-matrix |6| describing the one-dimensional Hubbard model Q. 

It soon became clear that for a full description of the string model, the R-matrices 
in higher dimensional representations corresponding to bound states were needed [s]. 
These bound state R-matrices could be computed directly by invoking the Yang-Baxter 
equation |9| or by using Yangian symmetry |10|. However, how to obtain these matrices 
directly from the fundamental R-matrix was unknown since the usual fusion procedure 
breaks down. In this note we will introduce a slight generalization of the fusion procedure 
for integrable models which allows us to obtain bound state R-matrices for the AdS/CFT 
S-matrix and Shastry’s R-matrix. 


This paper is organized as follows. In 

Section 2 

we 

introduce our fusion procedure and 

study it at the level of R-matrices. Then 

in 

Section 3 

we study more advanced aspects of 

it. Afterwards we turn to applications: In 

Section 4 

we discuss the well-known example of 


the XXX spin chain before moving on to the novel case of AdS/CFT worldsheet scattering 
and the one-dimensional Hubbard model in [Section ^ 


2 Fusion 

Consider an integrable system whose fundamental degrees of freedom are described by an 
n-dimensional vector space V^. These might represent the spin degrees of freedom of an 
integrable spin chain or the particle flavors of an integrable scattering problem in 1 -|- 1 


2 









dimensions. Their interactions are described by an x n^)-dimensional R-matrix 

7^(Ml,M2) : ^ ( 2 . 1 ) 

where the parameters mi,M 2 G describe the inhomogeneities of the spin sites or the 
particle rapidities. The space of parameters is typically a one-dimensional complex 
manifold]^ such as the complex plane or the Riemann sphere with certain punctures. 

For an integrable system, the R-matrix satisfies the Yang-Baxter equation and the 
involution property (the latter relation is understood up to an overall factor) 

7^.127^137^23 = 7^237^137^12, T ^ 12 T ^21 ~ li2- (2.2) 

Here and in the following, the indices denote the spaces in a tensor product as well as the 
associated parameters Uk- 

2.1 Singularities 

Suppose that there are pairs of points 

M(i2) = (mi,M 2) e M® C X M^, (2.3) 

where the R-matrix becomes non-invertible]^ In other words, the rank of the R-matrix 
drops below its maximum 

rank7?.(Mi, M 2 ) =(2.4) 

Commonly, these points form a one-dimensional sub-manifolc^ M® of x M^. The 
point M(i 2 ) € M® can thus be treated as a continuous parameter, and it is on the same 
footing as the parameters Uk G M^. In scattering theory, such singular behavior signals 
the existence of composite particles which are naturally part of the physical scattering 
problem]^ For these we will define new R-matrices to be interpreted as scattering matrices 
involving composite particles. We can thus extend the integrable system by adding these 
composite particles as additional degrees of freedom. 

More precisely, we are led to the introduction of an m-dimensional vector space V® 
by ‘fusing’ two spaces such that the extended integrable system is dehned on : = 
© V^. The corresponding R-matrix can be schematically written in block form 

^ext . yext ^ yext ^ yext ^ yext^ ^ext _ diag( 7 ^^^, 7 ^^^, 7 ^^^, 7 ^^^), (2.5) 

where the various blocks are maps of the types (H* = F, B) 

7^'^^^2(Ml,M2) : ^ (mi,M2) e X ( 2 . 6 ) 

^Additional continuous parameters or even discrete parameters are conceivable as well. 

^The R-matrix can have singular points where some of its eigenvalues diverge. As we are not interested 
in overall factors of the R-matrix (which may well depend on ui and U2) we should rescale the R-matrix 
at these points to remove the singularity. In other words, what counts is that the leading contribution in 
a Laurent expansion is non-invertible. 

^Discrete points with this property can also be considered along the same lines. We will, however, 
mostly be interested in the case of continuous values for U(i2) • 

^In fact, they are only part of the physical problem if they are also bound. The distinction between 
bound and unbound composite particles makes no difference here, and we shall refer to them collectively 
as composite states. 
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In particular, = TZ. The fact that satishes the Yang-Baxter equation is then 
equivalent to the statement 

^AiA2 ^AiAs^A2A3 'TPA 2 A 3 ^AiA^^A1A2 

/V12 /v^3 /V23 — /V23 /V^3 1-^12 

Moreover, the R-matrices have the involution property 

Of course the fusion procedure can be recursively applied, leading to larger and larger 
integrable systems. 


(2.7) 

( 2 . 8 ) 


2.2 Procedure 

In the following, we shall construct the fused R-matrices and afterwards check their prop¬ 
erties. 


Embedding and fusion matrices. We consider the R-matrix 77 (m(i2 )) := 7l{ui,U2) 
at a point M(i 2 ) where the rank drops to m < n^. At this point we can decompose it as a 
product of three matrices 

T(m(12>) : ^ V^®V^, 

'H(m(i2)):V^^V^, 

J^(m{12>) : ^ V®, (2.9) 


with the properties 


'^(m(i 2 >) = ^^(m(i 2 >)'H(m(i 2 >).T(m(i 2 )), J'(m(i 2 )) T(m(i 2 )) = 1®, 77 iuvertible. (2.10) 

Here, is a (surjective) m x matrix which fuses the tensor product ® to the 
space V®, and S is an (injective) n? x m matrix which embeds into ® V^. For 


convenience, we assume these two matrices to be pseudo-inverses as in (2.10). Finally, 77 
is a (bijective) m x m matrix and it represents the action of 77 on the space V®. In the 
following we shall describe two ways of obtaining the decomposition. 

First of all, let V® be the image of 77. Thus 77 can be understood as a map 0 
—)■ V®. Furthermore, dehne £ as the trivial embedding map V® —)■ ® V^. The 

above decomposition can be obtained as 7^ = (77T)“^77 and 77 = 77T. Note that this 
construction requires 77 = 77T to be invertible. This combination is not invertible precisely 
if the image of 77 contains a vector that is in the kernel of 77. In other words, the map 77 
contains a non-trivial nilpotent part. This case is more difficult to handle, and we shall 
exclude it for the time being. Later on in Sec. 5.3[ we shall discuss an explicit example. 

An alternative construction of the decomposition uses eigenvectors where we assume 
that the Jordan decomposition is trivial (i.e. the nilpotent case is excluded). The matrix 
77 possesses m non-zero eigenvalues Ai,..., Am- Let ci,..., Cm denote the associated right 
eigenvectors of 77 and e^,..., e™ the left eigenvectors 




e^n = Afee^ 


( 2 . 11 ) 


All of these quantities are functions of M(i 2 )- We normalize the vectors such that they 
form two dual bases for the space V® 

e^ei = 5l 77 = V A^e^e". (2.12) 
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Figure 1: Diagrammatical representation of fusion. 


We then define "H and T as the matrices of eigenvectors and eigenvalues 


E = (ei 


62 


e^) , "H = diag(Ai,..., Am), T = 


By their dehnition E and T satisfy the relations (2.10). 

The relations (2.10) imply the following useful identities 


^(M{i2))'H(n{i2)) = 7^(^t(l2))T(n(l2}), 
'^(M{ 12 )) -^(M{ 12 )) = -^(M{ 12 )) ^(M{ 12 >), 
EI{U(x2))E{U(i2)) -T('U(12>) = ^(M(12))- 


e^) (2.13) 


(2.14) 


These are the crucial relations that enable us to carry out the fusion procedure. 


Fused R-matrices. We introduce R-matrices by using E^ to ‘fuse’ together two spaces 
Vf 0 Vf into V^i 2 ) with Mp 2 ) = (wi,M 2 )- In particular, we are led to (cf. Fig. 1) 


"^fl2}3('^(12>, "Ws) -F(i2)(n(12)) Uz) 7^23(w2, Ms) ^(12)(M(12}), (2.15) 

"^™ 3 }('^l) '*^(23)) := ^ (23>(M{23)) ”^ 13 (^ 1 , M 3 ) ’^12(Mi, M 2 ) ^^{23} (M(23>)• (2.16) 


Any state from V^j^ 2 ) is niapped to Vf 0 V 2 ® V 3 by E^ acted upon with 7?. and then 
mapped back by T. 

Notice the different ordering of the R-matrices which is needed when spaces two and 
three are fused rather than spaces one and two. The R-matrix can then be dehned 
by applying the fusion procedure twice 


"^(12){34)('*^{12}, M(34)) -^(34) (M(34} ) 7^(12)4('*^{12) , M 4 ) 7^(12)3 ("^{IS} , M 3 ) i^{34) (M(34) ). (2.17) 


,bf 


BF 


In particular, from (2.14) it is readily seen that ( ]2.17[ ) can be cast into the symmetric 
form 


77-(i2}(34) ^{i2)^{zA)EZii{R-2iEZizEZ2zE{i2)E{zA)) (2.18) 

which demonstrates that it is independent of the order of fusing the underlying spaces. 
Here and in the following, we drop the parameters u and the labels F, B in favor of a more 
concise presentation. They can be recovered from the labels of the associated spaces. 

2.3 Relations 

In order to show that these R-matrices indeed describe an integrable system we have to 
show that they are invertible and that they satisfy the Yang-Baxter equation. 
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Involution property. We have to show the involution property (2.8) which reads more 
explicitly 


'^(12)3('*^(12})'*^3)’^3(12}('*^3;'^( 12 )) ~ 1 ; 


'^(12) (34) ) '^(34) )’^(34)(12) ('^{34}; '^(12)) 


BB 


1 . 


(2.19) 

( 2 . 20 ) 


Let us prove the hrst instance. For conciseness we will omit the arguments. Furthermore, 
we shall put brackets around the terms to be transformed in the next step 

^(12)’7^(12}3’^3{12) = [’H(l2)-Fp2)]7^13’^23^^(12)-^(12)’^32’^31^(12) 

= -^( 12 ) [’^ 12 ’^ 13 ’^ 23 ] £{ 12 )^ ( 12 )’^ 32 ’^ 31 <^( 12 ) 

= T (12)"7^23’^13 ["^ 12 ^ ( 12 )-^( 12 )] ’^32’^31^(12) 

= J-i^\2) [’^ 23 ’^ 13 ’^ 12 ]’^ 32 ’^ 31 ^^( 12 ) 

= \J-{I2)'ni2\ ["^13 [’^23’^32]’^3l]^(12) 

~ T-i{l2) [-^(12)<^^(12)] 

= H(i 2 ). ( 2 . 21 ) 

Invertibility of H then gives the desired result. In general, the strategy is to remove inter¬ 


mediate factors of E and T by use of the Yang-Baxter equation (2.2) and the properties 


(2.14) 


Yang—Baxter equation. Subsequently, it also follows directly from (2.14) that the 


above introduced R-matrices (2.15)-(2.17) satisfy the various versions of the Yang-Baxter 


equation outlined in (2.7) 


For example, since 77 satishes the Yang-Baxter equation (2.2) itself, (2.14) yields (we 
again suppress the explicit arguments) 

’^( 12 ) 77 ( 12 ) 377 ( 12)47734 = \H{12)^ ( 12 )] 'T^\‘{^2‘iE(^\2)^ ( 12 ) 77^7724 [^^( 12 ) 7734 ] 

= ^ ( 12 ) |]77i277i37723£^(l2)].F(i2)77i477247734£^(12) 

= J - (12)77 i 2 [77137723771477247734] £^(12) 

= [-7^(12)77127734] 77i4772477i37723£^(12) 

= 7734^^(12) [ 77 i 277 i 47724 ] 77 i 37723 £^( 12 ) 

= 7734 (12)7712] 77i47724£^(12)-F(12)77i37723£^(12) 

= [7734'R(12)] [-F(i2)77i47724£^(12)] [-F(12)77137723^^(12)] 

= 'H(i2 ) 773477(12)477(12)3. (2.22) 


This proves that the R-matrices (2.15)-(2.17) indeed describe the scattering of a new 


(composite) type of particle in this model. 


3 Further properties 

We have established the basic features of fused R-matrices. In the following we will discuss 
further properties. 
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3.1 Algebra 

Suppose there is a Hopf algebra H describing the symmetries of our integrable system. In 
particular, the R-matrix, by definition, intertwines the coproduct and opposite coproduct 
in the representation p^{u) : H —)• under which our fundamental degrees of freedom 
transform, i.e. for any generator JJ G H 


(pf(Ml) ® P2(“2)) = Tli2iUi,U2) (pf(Ml) 0^2 (“ 2 )) [A(J)] . (3.1) 

We dehne a new representation p^(m(i 2 )) : H —)■ V® for the composite degrees of freedom 
pfl2>(«{12>) [J] := .^(12>(M{12)) (pf (mi) ® P 2 (“ 2 )) [A(J)] ^(12>(M{12>)- (3.2) 


We will refer to this as the fused or composite representation. This representation clearly 
is m-dimensional. 

Let us show that (3.2) indeed dehnes a representation by proving that it respects 
the multiplicative structure. We have from (2.14) and the cocommutativity (3.1) of the 
coproduct, that for any J, J' G H 


’H(12)P(12) [J]P(12) [Jl — J^{12)T^12 {pi ® P2) [A(J)] £{12)J^{12) (pi ® P2) [A(J')] T(12) 

= ‘^( 12 ) (pi ® P 2 ) [A°'’(J)] 71 . 128 (12)^{12) (pi ® P 2 ) [A(J')] £{i2) 

= J^{ 12 ) (pi ® P 2 ) [A°P(J)] 7^l2 (pi ® P 2 ) [A(J')] 8 {i 2 ) 

= J^(12)’^12 (pi 8 ) P 2 ) [A(J)] (pi 0 P 2 ) [A(J')] T(12) 

= 'H(12).T(i 2) (pi 0 P 2 ) [A(J/)] T(i2) 

= 'H(12)P{12>[J/]• (3.3) 


Furthermore, the R-matrices (2.15)-(2.17) naturally intertwine the coproduct in the new 
representation. Explicitly, 


(p' 0 py) [A°P(J)] ipt' 0 P") [A(J)]. ( 3 . 4 ) 

For instance, let us prove the intertwining relation for the case Ai = B, A 2 = F. This is 
most conveniently done in the Sweedler notation A(J) = 0 J( 2 )- Co-associativity 

of the Hopf algebra is then written as 


^ J(l)(l) 0 J(l)( 2 ) 0 J( 2 ) = ^ J(l) 0 J( 2 )(l) 0 J( 2 )( 2 ) 


and the intertwining property (3.1) of the R-matrix is formulated as 
^77.12Pi[J( 1)] P2[J(2)] = y^Pl[J(2)] P2[J(1)] 7li2. 


(3.5) 


(3.6) 


In this language, we have (for conciseness we suppress the arguments u, v of the R-matrices 
and the sums) 


’^(12)3(p(12) 0 Pa) [A(J)] = 77 .(i2)3J^(12) Pi[J( 1 )] P 2 [J( 2 )( 1 )] P3[J( 2 )( 2 )]^^( 12 ) 

= 7F (12)77.1377.23 Pl[j(l)] P2[j(2)(l)] P3[j(2)(2)]^^(12) 

= .F(i2) Pi[J( 2)(1)] P2[J(2)(2)] P3[J(l)]^^(12).F(i2)77i37723i^(12) 

= (p(i 2 ) 0 Pa) [A°p(J)] 77(12)3, (3.7) 


where we used (2.14) repeatedly and (3.6) in the third step. This proves that 77^^ displays 
the expected cocommutativity properties. 
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3.2 Auxiliary features 

Here we will discuss some auxiliary features of the fused R-matrices. 


Similarity transformations. We have the freedom to apply a similarity transforma¬ 
tion W{u(i 2 )) to the space V® 

^•('^( 12 )) £^(M( 12 )) 1 R(U( 12 )) 

^{U( 12 )) —W {u ( 12 )) J-{u ( 12 )) 1 

'H(M(12)) -t w {U(^12))U{U(^12))W {U(^12))~^. (3.8) 


This transformation affects none of the relations (2.10), and therefore all the above results 
apply to the transformed system without further ado. In the construction of 8,71, T via 
eigenvectors (2.13), the similarity transformation does not preserve the diagonal nature of 
Ti (unless W is diagonal as well). However, we have not made explicit use of this property 
in the constructions. 


Opposite form. A similarity transformation by 'H(i 2 ) has a curious effect on the fused 
R-matrices 


'^(12>^(12)3'^(’i2) = ’^(12).^(12)’^13^23^^(12)'^(12) 

= vT(12)’7^12’^13’^23^(12)’Hp2) 
= { 12 ) 7 ^ 237 ^ 137^128 { 12)71 ^^ 2 ) 
= J^{12)7l237ll38{12) ■ 


(3.9) 


Compared to the original dehnition 7^(i2)3 = 7 ^{i 2 ) 7 ^i 37 ^ 238 {i 2 )-i we observe that conjuga¬ 
tion by 7 i{i 2 ) flips the order of the R-matrix factors within the fused R-matrix. 

This observation goes hand in hand with the definition (3.2) of the fused representation 


via the coproduct. If instead of the coproduct we use the opposite coproduct, the resulting 
representation is related to the original one by a simple similarity transformation 


J^{ 12 ) [pl ® P 2 ) [A(J)] 8 { 12 ) - 7 ii^^ 2 ) -^(12) (pi ® P 2 ) [A°P(J)] 8 { 12 ) 7 i{l 2 ). 
As usual we used (|2.14[) and the intertwining property of the R-matrix. 


(3.10) 


Symmetric R-matrices. In many practical applications 7 ^(m(i2 )), acting as an operator 
on the space ® V^, is symmetric w.r.t. some inner product, e.g. the standard inner 
product (a, 6) := aJh defined on V^. In most cases, 7 Z{u{i 2 )) = 7 ^(m(i2 ))''' admits an 
orthonormal basis of eigenvectors]^ In addition, the fusion and embedding matrices are 
conjugate to each other T = 8 ^. 


One minor problem is that the resulting R-matrices (2.15) are not symmetric. Trans¬ 


position reverses the order of the constituent R-matrices which corresponds to a similarity 


transformation according to (3.9) 


’^(12)3 


= .^( 12 )’^ 23 ’^ 13 ^^( 12 ) = 7i(i2\7l(i2\37i 


( 12 ) /<-(12)3^12^ • 


(3.11) 


®This is evident if the inner product is positive definite. For indefinite inner products (including 
complex symmetric matrices), eigenvectors can be null. In this case the eigenvectors cannot be normalized, 
and even a non-trivial Jordan decomposition may arise. 










By applying a similarity transformation (3.8) defined by a ly such that "H = W'^^W, the 
resulting R-matrices become symmetric 


Km = = w-Ki2}'^(mK2y^ = y-^Km^K = Km 

(3.12) 


3.3 Complementary fusion 

A fused R-matrix can also be dehned for the complement of the space V^i 2 ) 

Vf ® = vf,2) © vf^^, Km ■ ^2) ® ^3 ^ ® (3-13) 

As we shall see, this space is even better suited for fusion. 


Complementary space. To dehne an R-matrix for the complement, we need to dehne 
fusion and embedding matrices for the complementary space 

T(m(i2)) : ^ © V^, ^(m{12>) : ^ V^. (3.14) 

They are dehned to obey the following orthogonality properties with the original embed¬ 
ding and fusion matrices in ( 3.15| : 

J'(M(12)) £^(M(12>) = 1^, J^{uii2))8{uii2)) = J'(M(12}) (M(12>) = 0. (3.15) 

They directly imply the completeness relations 

^(12).y(12> + ^^(T2>-^{T2> = ll2- (3.16) 

as well as orthogonality relations with the R-matrix 

Kiu{12)) S{u{12)) = 0, {u{12)) Tliu{12)) = 0. (3.17) 


Complementary R-matrix. One can dehne a complementary R-matrix in analogy to 

'^^>3('^(12>) "^s) '■= ^ (T2)('^{12)) 77.13(Mi, M 3 ) 77.23(-U2, %) ^^{T2)('^{12))- (3.18) 

The other related R-matrices follow in a similar fashion. The various integrability rela¬ 
tionships can be derived in a similar fashion to the above. The general strategy is to 
remove the factors of T^X 2 }'^{T 2 } which typically appear between the various R-matrices. 
The starting point is the relationship {i 2 )'nu'R- 2 z£{^) = 0 which follows from the above 
dehnitions 


'H{i2)J^{l2)T^l^'T^2‘i£{12) — .^( 12 ) 7 ^ 127 ^ 13 ’^ 23 <^( 12 ) ~ .^( 12 ) 77 . 23 ’^ 13 ’^ 12 ^^( 12 ) ” 0 - 


(3.19) 


Together with the completeness relationship (3.16) one can show 


£{12)^ {12)"^13’^23<^(12) ” ’^13’^23^^(12) ~ ^(12).^(12)7^13’^23^^(12) ” ’^13’^23<^(12)• (3.20) 

Therefore, all intermediate factors of £{^)J^{\ 2 ) can be removed iteratively from products 
of R-matrices from the left to the right. 
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Similarly, one can prove that 


^ ( T 2>^2 


-'23’^13^(12} —0, ^^(12)-^(12)’^23’^13^(12) — ^23^13^(12)- (3-21) 

On the level of the R-matrix this simply corresponds to a similarity transformation on 


the composite particle space applied to the original fused R-matrix (3.9). 


Put differently, the combinations 7^137^23 and 7^237^13, when viewed as a block diagonal 
matrix, effectively possess a triangular shape: 




13 'V23 


>fi2)®Vf ^(Vf,2)©V|5))®Vf, 


{ 12 > 


( 12 > 


7^237?. 


23 'Vi3 


V|2) © Vf ^ Vfi 2^ 0 


(3.22) 


Composite states from the complementary subspace are mapped to the same sub¬ 
space by 7^.137^23• Conversely, composite states from the original subspace V® 2 ) can map 
to both subspaces. Fusion as dehned in (2.15) works only due to the presence of the 
projector -F(i 2 }- However, for the opposite R-matrix (3.9) this reverses and states from 
V^j^ 2 ) are mapped to the same subspace. This is at the cost of a similarity transformation 
on Vf, 2 ). 


Opposite fusion. The involution property (2.2) also has an interesting implication on 
the complementary space V® which we shall discuss in the following. According to our 
assumptions, the hrst m eigenvalues of TZ{u(^i 2 )) are non-zero while the other — m 
vanish at the point M(i 2 ) = {ui,U 2 )- We assume that the latter fall off linearly in e when 
approaching the singular point as in {ui,U 2 -|- e) —)■ M(i 2 )- Furthermore, also the product 
7^127^21 is assumed to be proportional to e. 


By the involution property (2.2) we have that away from the singular points 7 li2 
is invertible with inverse proportional to 7^21- Therefore 7^i2 and 7^21 share the same 
eigenvectors, but with inverse eigenvalues. Now, by our assumptions on the behavior near 
the singular point, we are led to the conclusion that the hrst m eigenvalues of 7^21 fall off 
linearly with m while the remaining v? — m eigenvalues approach a constant^ 

Therefore, the null space of 7^21 (w( 2 i)) at M( 2 i) := (n 2 ,Mi) is given by the image of 
^( 12 ) (^( 12 )) whereas the non-trivial eigenspace is the image of £^( 2 i)(m( 2 i})- Consequently, 
the space V® is the complement of V^i 2 } within Vf 0 Vf, and it has dimension ffi : = 
— m. We therefore reproduce the elements of complementary fusion 


vh , — vh_ 
^(21) — 








(T2)’ ~ -^(12>) ^(21>~^{12>- (3.23) 

All the constructions for the complementary fused states and operators proceed as before 
with the roles of spaces 1 and 2 interchanged. 


3.4 Algebraic Bethe ansatz 

Let us briehy touch upon the effect of fusion on monodromy and transfer matrices that 
play a central role in the algebraic Bethe Ansatz. 


®If some of the eigenvalues of 7?.i2 scale with a higher power of e (consequently also Tii 2 'R- 2 i), only 
fewer than r£ — m eigenvalues will be finite. The following considerations would have to be adapted 


accordingly. See Sec. 5.3 for a concrete example of this case. 
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RTT-relation. The key relation in the algebraic Bethe ansatz is the so-called RTT- 
relation. The RTT-relation describes the commutation relations between the elements of 
an n X n dimensional operator valued matrix T^('u) : —)■ called the monodromy 

matrix. We have 


nll{u^,U2)T^{ur)Ti{u2) = Ti{u2)Tl{u{)nll{u^,U2). (3.24) 

However, rather than taking 7?.fJ(ni,n2) one can also consider taking fused R-matrices 
and consider the RTT relation this would generate. To this end, we introduce a fused 
monodromy matrix 

^?2>(“(12>) '■= ^{12){U(^2 ))Ti {Ui)T2 {U 2 ) £{12){U(^2))- (3.25) 

It is straightforward to show that 'T^ 2 ) satishes the RTT-relation for fused R-matrices. 

An object of special interest is the transfer matrix, which is dehned as the trace of the 
monodromy matrix 

t^iu) = trp T^(n), ^^(w{i 2 >) = Rb T®(m(i 2 >)- (3.26) 

This object generates the mutually commuting set of operators that is the dehning prop¬ 
erty of integrable systems. 


Fusion for transfer matrices. We can formulate the relation between the transfer 
matrices in different representations. Consider the product of two transfer matrices and 
use the completeness relation ( |3.16[ ) 

f{ui)t^{u2) = tri2[7j(Ml)7^(M2)] 

= trp 2 ) [-^(12)’7i(Mi)72(m 2)^(12)] + tr( 21 ) \_^{ 21 )'Tl{Ui)T 2 {u 2 )S{ 2 l)\ 

= t''(w(i2>)+tV{2i>). (3.27) 


For the latter term, we note that a similarity transformation by 77(21) is required to 
interchange the order of monodromy matrices in analogy to (3.9). 


4 The Heisenberg XXX spin chain 


To illustrate our fusion procedure, let us hrst consider the Heisenberg XXX spin chain. 
The fundamental particles transform in the spin-| representation of su(2). The corre¬ 
sponding R-matrix is of difference form 77xxx('Ui, ^ 2 ) = 77xxx('Ui — M 2 ) and is given by 


77xxx(m) 



+ 1 

0 

0 

0 

\ 


0 

u 

1 

0 



0 

1 

u 

0 


\ 

0 

0 

0 

u + 

1/ 


(4.1) 


There are two points at which the rank of TZxxx{u) is not maximal, namely at u = ±1. 
At M = 1 the rank is three, while at w = — 1 the rank is one. These points correspond to 
the singlet and triplet that arise in the decomposition of the tensor product of two spin- 7 
represent at ions. 
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Singlet. At the point ui — U 2 = —1, there is only one eigenvector with non-zero eigen¬ 
value and our projector is 


E = r = — 


1 


/0\ 

1 

-1 

VO/ 


U = - 2 . 


(4.2) 


It is easy to check that the identities from (2.14) hold. The representation that is gen¬ 
erated according to (3.2) is the trivial representation. Consequently, the fused R-matrix 


describing the scattering of a singlet state with a fundamental particle is given by 


’^(12)3(M(12>, Ms) — [{U2 “ Ms)^ “ l] 


1 0 
0 1 


(4.3) 


In other words, we see that the singlet has trivial scattering with a doublet (up to an 
overall factor which depends on the dehnition of the overall factor of T^xxx)- 


Triplet. The other point ui — U 2 = +1 is the opposite of ui — U 2 = —1. The resulting 
space therefore is the complement of the above singlet. More concretely, there is only one 
null vector and from the three remaining eigenvectors we hnd 


S = J^^ 


/I 

0 

0 

VO 


0 0 \ 

0 1 / 


= diag(2,2,2). (4.4) 


It is easily seen that this will give rise to the usual spin-1 representation of 5 u( 2 ). Indeed, 
(|3.2[) yields for the simple roots 


/O 72 0 \ 

P{12>[S1= 0 0 72 , 

Vo 0 0 / 


/ 0 0 0 \ 

P(i2)[§-]= 72 0 0 . (4.5) 

V 0 72 0/ 


Consequently, we recover the standard 
triplet-doublet case from (2.15) 


composite state R-matrix 7 ^(i 2 } 3 (m(i 2 },M 3 ) for the 



/M2 

— M3 -|- 2 

0 

0 

0 

0 




0 

M2 - M3 

72 

0 

0 

0 

(mi - M3) 


0 

72 

M2 — M3 -l- 1 

0 

0 

0 


0 

0 

0 

M2 - M3 1 

72 

0 



0 

0 

0 

72 

M2 - M3 

0 


V 

0 

0 

0 

0 

0 

M2 - M3 + 2/ 


(4.6) 


One can now easily check that the Yang-Baxter equation holds. 


5 The Hubbard model 

The R-matrix for the Hubbard model is best described in terms of its symmetry algebra; 
centrally extended su(2|2). This algebra is obtained from regular su(2|2) by adjoining 
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two additional central charges to it. The R-matrix in the fundamental representation 
is completely hxed by the intertwining property (3.1) 011|. The R-matrices involving 


composite state representations are hxed by Yangian invariance 110 


It turns out that there are two cases where the R-matrix becomes of lower rank. There 
are two cases where it becomes of rank 8 and the fused representation is a (a)symmetric 
short representation. Moreover, there is a point where the R-matrix reduces to rank 1 
corresponding to a singlet representation. 

We will show that the bound state R-matrices found in the literature follow from 
our fusion procedure. 


5.1 The Hubbard model R-matrix 

The R-matrix for the Hubbard model is a 4^ x 4^ matrix. It acts on the tensor product of 
two 2|2-dimensionaI spaces with bosonic basis vectors |0“), a = 1,2 and their fermionic 
counterparts a = 3,4. For convenience let us introduce the 212-dimensional basis 

vector = (| 0 ^), | 0 ^), ^md let E^b be the matrix unities with a (— 1 ) 1 '®! in 

row A and column B. The fundamental R-matrix is then of the form 

TZ{uuU 2 ) = (-1)I^I+I^Ie^b 0 E^dTI^a^c{ui, ua) (5.1) 

with the only non-zero entries given by 


/yt I /yt I rf> ^y» /y» ”| 

a./1 a./0 vO 1 Uy-I JyQ -L 


'Tya c rare I /rare ro-r^N*^! *^2 '^1 *^1*^2 

n b d = 0^0,, +{df,da-0^0^)— -, 

/y ‘ /y ‘ /y /y I 

^ •jO‘2 vO 2 \ 


= 


U 2 xf 


Xr, 


Ui — X 2 




ry» I ry* ry* ry» 

*^2 *^2 1 


ry> I ry» /y ‘ /y ry> 

vt/1 '^‘2 ^^2 ^^2 


- 1 


a ^ ^ ^aj3 _j_ ^ _ 

a b ^ab^ 


I 1 I 2 U 2 X 1 X 2 {xt -xt) 


X, 


ry>~^ /y~' ( ,y /y T ^ /y~^ 

Xi X 2 ^2 ^^^2 

/ ry>~^ /y~^ / ry> / /y >^”1” 

I3 \^l ^2 )\^2 ^2) 

7172 (^ 1 ( 11^2 - i)(^t - A) 


(5.2) 


and 


06 rarck 


1/y~^ ry>~^ 

vt/1 U/n 


Ui Xi 


X 


+ ’ 


7^“ 


y raro: '^2 
b = — 


Ih X 2 -x^ Ji 

Ui X 2 - (c/ 72 


7^ 


/v» I 

06 a raro^^^l 

b 0 - Obd/s^ 

Jbr\ 


^ 72 

r 7i 


X 


= 5t5^pU2 ^ 


Xr, 


Xn 


(5.3) 


The parameters xf 2 sxe related to the spectral parameters ui^ 2 , respectively in the fol¬ 
lowing way 

+ I h _ 1 h , 

u = x^ + — - - = x -I-h-. (5.4) 

x+ 2 x~ 2 

The parameter h corresponds to the coupling constant and the parameters U are related 
to the above by = x^/x~. Finally, the additional parameter 7 dehnes the relative 
normalization of bosons (0) and fermions ("0). The R-matrix is the intertwiner of the 
centrally extended 5 u( 2 | 2 ) superalgebra. This algebra contains two 5 u( 2 ) subalgebras 
that act on the bosons and fermions respectively. The R-matrix is a symmetric matrix 
with respect to an appropriately chosen inner product for the states 


(riy) = s 


ab 


u 


( 3 .+ _ 3 .-) _ (JJ _ 


7 


T 


(5.5) 
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5.2 Symmetric states 

Let us first consider the points where the R-matrix becomes of rank 8 . They correspond 
to the two special points = X 2 (corresponding to Ui = U 2 + h) as well as x^ = X 2 
(corresponding to Ui = U 2 — h). In the following we will consider the point xf = x 2 ■ The 
considerations for the other point are completely analogous. 

This hrst ingre dien t we need is the matrix of normalized eigenvectors £ of the funda¬ 
mental R-matrix (5.1). There are four bosonic vectors and Three of them 

have unit eigenvalue \ab) = 1 and form a standard triplet under the bosonic 5u(2) 




= ir> ® ir), 

and the fourth contains two fermionic states and is a singlet under both su( 2 )’s 
|g[34|, ^ 7l72t/2'‘ ® \i’^)] + (1 - W'l" + Uj)) 

- xr) + (1 - 1 X 2 "- + uT}f 

rf> 

vO 1 vL 2 


( 6 . 6 ) 


-^[34] — 


X 2 X 2 

1 — XiXi 


'2 \-^2 
^2 ~ 


+ 


^2 


UiU, 


2 Xn 


X, 


There are four fermionic eigenstates |S““) 

72 |<^“) <8) 1^/’") 7l^2 |^“) <8) |0“) 


|S““) = 




1 /y> 

vL 9 R/ 9 


Xn 


+ f/l4 


Xi Xi 


Xi 


Xn 


(5.7) 


(5.8) 


\/U2{x2 - X 2 ) + UiU^ixf - X^ ) ’ U 2 Xi 

After packaging these vectors in the matrix £ and taking = T"'", it is readily checked 
that (2.10) and (2.14) hold. 

We are now in a position to apply our fusion procedure. Let us consider the fused 
matrix (2.15) and make it symmetric by applying the transformation W as in ( 3.12[ ). It 
is then straight-forward to check that this R-matrix coincides with from |[^|. 

As a non-trivial example let us compute two scattering processes 


77'^^(M(12>,M3)|S^^^^) ® 10 ^), 




13\ 


( 5 . 9 ) 


It is easy to see that 

77'^"^(m{i2},M3)|S^^^^) ® 10^) = ( 8 ) 10^), (5.10) 

which agrees with oi = 1 from section 6.1.2 of |9|. This shows that we have the same 
normalization for the bound state S-matrix. Subsequently, we have by dehnition 


® 10 ^) 

= WJ-i277i37723^12lR"^|S'=') ® 10=*) 

^ ' U 2 ixt - X^){xt - 3^3 ) ( 3^3 - ^1 
(x^-x^)(x)|'-X 2 ) 

^ U 2 {.xt - X 2 ){X^ - X-t){x-^ 

( 'T* '' 

yxi -^3 )\-^3 -^ 2 ) 

U 3 


U 2 {xt - X 2 

U3 xt - X3 

U 1 U 2 x^ — x^ 


) + t/iR|(x+ 


-Xi 


_|_ \^2 ^2 )\^3 '^1 )\'^3 ^ 2 ) 

Ui{x^ -xt){x2 -x^) 

- 4 )' 


|S^=*) 0 10 =*) 


(5.11) 
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where we used the inner product ( |5.5[) and = X 2 . This result exactly coincides with 
the literature, in particular it is the coefficient ay from section 6.1.2 of P|. In fact it is 


not hard to check that the representation (3.2) is exactly the two particle bound state 
representation from ||^|. 

Let us conclude this section by considering the complementary fusion procedure. The 
complement is spanned by antisymmetric states |A) = |A““)} that are 

perpendicular to the symmetric states |S) given above. 


For these states, one can then easily compute the complementary R-matrix (3.18) 


In particular, due to the upper triangular structure ( 3.22[ ), the computation simplihes 
somewhat. For instance, let us consider the analogue of the states (|5.9[) 


|A( 33 )) (g, 1 ^ 3 ^^ 


|A^3) (g) 10^). 


(5.12) 


The action of 7 ^ 137^23 on both states is simply multiplicative and yields 


xT - X, U 3 


■Rl3K23|Al”>) 0 ^ 

Xi — xj Ui 


X, U3 


^2 -3^3 U 2 


|A(33)) 0 


xj - X, U3 


7^137^23|A'") ® 10') = 4- 

Xi — xj Ui 


Xn 


X3 U 3 U 1 U 2 Xi 


X, 


-x+ U 2 


U 3 


(2/0 Xq 


(5.13) 

lr|Ay®|y). (5.14) 


This then corresponds to the fused R-matrix on antisymmetric states corresponding to the 
other point of lower rank x^ = x^. Notice that the matrix has a different normalization 


and the inverse of the coefficient (5.11) appears, which is in agreement with |12| where 


the relation between the symmetric and antisymmetric R-matrices is discussed. 


5.3 Singlet state 


The R-matrix coefficients in (5.2) have a common factor of (x0 — 4 )/(^r ^2 ~ which 
has a potential singularity at x^ = l/x^j^ Furthermore, x^ 7 ^ x0 in order for the 
numerator to be non-zero, i.e. x^ = l/x),*". The singularity affects only the su(2) x 5 u( 2 ) 
singlet sector spanned by the two states 


|BB) = k 


ab , 


IFF) = £“'^ 4 '“ 0 V'* 


(6.16) 


The action on the remaining 14 states is hnite. Acting on the singlet, the R-matrix reduces 
to a 2 X 2 matrix M. In the limit xf —)■ l/x)*", it has the following singularity structure 

M = + . (5.16) 

Up to an overall factor we hnd for the residue 




(uif/27i72/(f/i - f/f') 


-(Ui - t/r')/7i72\ 

-U 1 U 2 J 


(5.17) 


We know that = x^/x~ and hence U 1 U 2 = ±1. Both values of the latter signs are 
permitted, and we have to discuss the two cases separately, as they lead to rather distinct 
behavior. 

^Here, we do not explicitly multiply the R-matrix by an overall factor to compensate the singularity. 
Hence we will consider the most singular contributions to the matrix. 
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We start with U 1 U 2 = —1 which is analogous to the cases discussed above. There are 
two eigenvectors 

(t/i - t/r')|BB) + 7172 IFF) and {U, - U{^)\BB) - 7172 IFF). (5.18) 

The eigenvalues are 2 and 0, respectively, therefore only the former state is singular. It is 
a singlet of the Yangian algebra, however it has a non-trivial charge lF(i 2 ) = U 1 U 2 = —1 
which plays a role for the coproduct. The other state belongs to an adjoint representation 
of psu(2|2)|^ The singularity structure of this case is peculiar with respect to comple¬ 
mentary fusion: In the limit xf —)■ l/xj*", the eigenvalue of the above non-singular singlet 
state approaches zero even faster than for the 14 non-singlet states. We hnd 


Ala 


A 


14 


Alb ~ e. 


(5.19) 


This means that complementary fusion based on 7^21 produces merely one composite state 
rather than 15. This state is just the other singlet, and in fact one can see that exchanging 
the two sites in ( 5.18[ ) interchanges the two states|^ 

The other case 17(12} = U 1 U 2 = -l-l has only one eigenvector 

(17i-f/fi)|BB)+7i72|FF), (5.20) 

and its eigenvalue is 0. The matrix does not admit a second eigenvector because it 

has a non-trivial Jordan decomposition. This makes the case very special, and the fusion 
procedure described in this paper does not immediately apply. Let us therefore try to 
understand what is going on: The eigenstate is the singlet state discussed in |4|. There 
it was shown that the state behaves just like a fused state under scattering with other 
states, i.e. it preserves its form. To understand the role of the other singlet state, it makes 
sense to take a closer look at the limit xf —)■ l/x^. Here, both singlet eigenvalues remain 
finite, whereas the eigenvectors become collinear. Therefore, at xf = 1/xf there is only 
one meaningful eigenvector, the difference of the eigenvectors plays no signihcant role. 
Even though all eigenvalues remain hnite at this point, fusion does take place due to the 
coincidence of eigenvectors. 

The fused R-matrix 7^(i2)3 of the singlet state with 17(i2) = -|-1 was shown to be trivial 
Q up to an overall phase factor related to crossing symmetry 


m 
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Based on this result 

one can easily derive the fused R-matrix of the singlet state with f/(i 2 ) = — 1- The point 


is that the factor Ui appears in odd powers in (5.2) and (5.3) only when the third index 


is fermionic. Flipping the sign U 2 flips the sign of 7 ?.(i 2)3 precisely if state 3 is fermionic 
( |5.1[ ). This means that the state with 17{i2) = — 1 behaves like a fermionic singlet, while 
U{i 2 ) = +1 corresponds to a bosonic singlet. This observation is in line with the coproduct 
rule of odd generators. 
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